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A considerable understanding of the formal description of quantum mechanics has been achieved 
after Berry's discovery |2| of a geometric feature related to the motion of a quantum system. 
He showed that the wave function of a quantum object retains a memory of its evolution in its 
complex phase argument, which, apart from the usual dynamical contribution, only depends 
on the "geometry" of the path traversed by the system. Known as the geometric phase factor, 
this contribution originates from the very heart of the structure of quantum mechanics. 

A renewed interest in geometric phenomena in quantum physics has been recently mo- 
tivated by the proposal of using geometric phases for quantum computation. Geometric (or 
'Berry') phases depend only on the geometry of the path executed, and are therefore resilient 
to certain types of errors. The idea is to exploit this inherent robustness provided by the 
topological properties of some quantum systems as a means of constructing built-in fault tol- 
erant quantum logic gates. Various strategies have been proposed to reach this goal, some of 
them making use of purely geometric evolutions, i.e. non-Abelian holonomies |29 20| |19I . 
Others make use of hybrid strategies that combine together geometrical and dynamical evo- 
lutions 1111 . and others yet use more topological structures to design quantum memo- 
ries (5][21. Several proposals for geometric quantum computations have been suggested and 
realised in different contexts, in NMR experiments 1 11 1, ion traps f6l ll4ll26lllOI 1271 . cavity 
QED experiments |23|, atomic ensembles |28, 15 1, Josephson junction devices |8|, anyonic 
systems |13|, quantum dots 1231. 



Geometric phase and holonomy 

Suppose that a system undergoing a cyclic evolution is described by classical mechanics; it is 
impossible to tell from its initial and final state whether it has undergone any physical motion. 
The situation in quantum mechanics is quite different. The state vector of a quantum system 
retains the "history" of its evolution in the form of a geometric phase factor. 

This deep and fundamental concept was originally discovered by Pancharatnam I 2D in 
the context of a classical beam of polarised light and "rediscovered" in a quantum mechan- 
ical context by Berry |2|. Pancharatnam introduced the concept of parallelism between two 
states, as a criterion to compare the relative phase between two beams of light with different 
polarisation. He recognised that a natural convention to measure the phase difference be- 
tween two interfering beams is to choose a reference where the intensity has its maximum. 
For example, by superimposing two beams of polarisations ipi and ip2 the intensity is propor- 



2 



tional to / oc 1 + \ \tp2)\ cos (x + arg (tpi \ip2))- The interference fringes are shifted by 
(f = arg {tpi \ip2)i which, following Pancharatnam's prescription, represents the phase differ- 
ence between ipi and -02. This idea, translated into quantum mechanics, leads to the definition 
of relative phase between any (non-orthogonal) states lying in a (finite or infinite) Hilbert 
space. When arg (tpi \ip 2 ) = 0, tpi and ip2 are called in phase. Pancharatnam's most impor- 
tant contribution was to point out that this condition is not transitive: if tpi is in phase with ip2 
and ip2 with ips, the phase between ipi and ip^ is, in general, not zero. As in quantum mechan- 
ics states are defined up to a phase, tp 2 can always be redefined parallel to ipi . However, when 
a third state ip3 is considered, it is, in general, impossible to redefine it in phase with both ipi 
and ip2- This is due to an irreducible phase contribution \ = arg (ip±\ ^2) (^l ^3) {^3 1 ^1), 
called Pancharatnam phasePancharatnam phase, which represents the most elementary exam- 
ple of geometric phase fTl fT7ll22l . 

If, instead of a discrete collection, we consider a continuous chain of states | (j)(s)) (with 
s E {so . . . si}), we can repeat a similar argument and redefine the local phases | 4>(s)) — > 
I if)(s)) = e la ^ \ 4>{s)) to impose the phase condition between infinitely neighboring states, 
namely 

arg (iP(s)\ ^{s + da)) ~ (ip(s)\ ^(s))ds = 0, (1) 

as 

which is known as parallel transport condition. As emphasized earlier, this condition is not 
transitive. Therefore, although neighboring states are in phase, states far apart along the curve 
accumulates a finite phase difference between them. In particular, if the chain is a closed loop, 
i.e. I (f>(s )) = I 0(si)), a state "parallel-transported" around the loop experiences a phase 
shift 

|V(»i)> =e < *M« )>, X^=a{si)-a{s 2 ) = iJ S \ct>\^- s \ct>)ds = ij> {<j>\d<j>), (2) 

which is the celebrated geometric phase. As for the Pancharatnam phase, x-y is an irreducible 
phase contribution which solely depends on the closed path 7 traced out by | tp(s)) in the 
Hibert space. It is easy to verify that neither a local redefinition of phase, nor a change in the 
rate of traversal affects the value x-y- 

Adiabatic implementation of holonomies 

A natural question to ask is how the idea of parallel transport applies to physical scenarios. It 
turns out that this concept plays a key role in a variety of physical contexts (see l24l l3l lT8l '). 
and, in quantum mechanics it emerges as a natural feature of adiabatically evolving systems. 

Suppose that an Hamiltonian, H(X t ) is controlled by a set of time-dependent parameters 
At. If the requirements for the adiabatic approximation (see 116111 21 ) are satisfied, a state, ini- 
tially prepared in an eigenstate | 4>(to)) = I ^n(M ))> remains eigenstate of the instantaneous 
Hamiltonian, during the evolution: 

\^(t)) =e 45 W|*„(A t )) where H (A 4 )| ¥ n (A t )> = e„(A t )| f n (A t )>, (3) 

(h = 0) where S(t) — — J t e n (X t )dt is the usual dynamical phase. Under this approximation, 
the state | ip{t)) can satisfy the Schrodinger equation only if the contraint {^ n {X t )\ ^^(Af)) = 
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is fufilled. Hence, the state | 5 , „(A t )) is parallel transported around the Hilbert space as the 
parameters A's are varied. If the latter are eventually brought back to their initial values Ao, 
and the eigenspace of | is non-degenerate, the final state will be proportional to the initial 
one, | ^(tf)) = e lx ~« \ ^(ij)), with an accumulated geometric phase x-y (which in this context 
is called Berry phase), only dependent on the path, 7, traced in the parameter space: 



where the path integral here is explicitly expressed in terms of a vector (one-form), known 
as Berry connection. The inherently geometric nature is even more evident when the path 
integral in eq. (0} is formulated as a surface integral, via the Stokes theorem: 



where £ is the surface enclosed within the loop traced by A in the parameters' manifold, and 
Fij — diAj — djAi is called the Berry curvature. The Berry curvature in many interesting 
cases (such as for qubits) is a slowly varying function, or even a constant. As a result of 
this, the geometric phase behaves as an area and depends almost exclusively on the surface 
enclosed by the loop. This is one of the crucial characteristic that makes the geometric phase 
quite appealing for the implementations of fault-tolerant quantum computation. A feature, 
like an area, which is much less dependent on the details of the time evolution, is likely to be 
less affected by variations of environmental conditions, and hence, more robust. 

The prototypical example in which this area-like behaviour is manifest, is the case of 
a single qubit adiabatically evolving under a generic Hamiltonian H{t) = n t • a, where 
<f = (cx,&y,&z) is the vector of Pauli matrices, and ri t is a time-dependent vector. It is 
possible to show that the curvature associated with a qubit state gives rise to a very simple 
form of the geometric phase, namely x-y = i§ (i depending on whether the qubit is initially 
aligned or against the direction of n), where is the solid angle spanned by the direction of 
the vector n. The curvature in this case is constant (±1/2) and is the surface enclosed in 
parameter manifold (the Bloch sphere) ( see Fig.^). 

Before turning the discussion towards the implementation of quantum computation, it is 
important to introduce the non-Abelian generalisation of the geometric phase, or holonomy. 
In obtaining the geometric phase for an adiabatic evolving system, the assumption that the 
eigenspace to which the prepared state belongs is non-degenerate was crucial. Such a con- 
dition insures that, when a loop in the parameter space is traversed, final and initial states 
are proportional: i.e. the net effect of the evolution is merely a phase. However, assuming a 
degenerate eigenspace, opens up a wider variety of possible evolutions, with a slightly more 
complex structure, known formally as holonomy. 

The word holonomy refers to the set of all the closed curves, or loops on a manifold, start- 
ing and ending in the same point xq . It is easy to verify that this set has the structure of group 1 . 
The geometric phases themselves form a representation of an holonomy group: any loop in the 

'The composition of two loops and is obtained by joining the end point of one loop with the starting point of the 
other. The identity element is the trivial loop with only one point (xq). The inverse of curve is the same traversed in 
the opposite direction. For a rigorous definition see II Sll9l ) 
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(5) 
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parameter space of an Hamiltonian is associated with a geometric phase factor. And clearly 
they form an Abelian representation as phases commute: e* XT i e 1 *' 1 ? = e* x *" e !X ^i . This 
therefore implies that their non- Abelian generalisation are not represented by ordinary num- 
bers, but by matrices. This naturaly emerges in adiabatic evolving systems, when eigenspaces 
are degenerate. 

Let's write a parameter dependent Hamiltonian in the form: H(X t ) = J2 k € k{Xt)^k{Xt), 
where ITfc(At) are the projector operators of the instantaneous eigenspaces. As time varies, 
the parameters change and with them eigenvalues and eigenspaces. The latter are smoothly 
concatenated via a unitary transformation O(Xt) (the eigenspaces never change dimension, 
as this is forbidden by the adiabatic requirements), II/.(At) = 0(X t )H k O' (Xt), where n" is 
an eigenspace at the initial time to (O(Xt ) = 1). The unitary transformation produces 
the change of picture to the frame moving rigidly with the instantaneous eigenspaces. In this 
frame, the evolution is governed by the Hamiltonian H(Xt) = ^2 k €hJl k — i d °dt O^i^t)- 
Imposing the adiabatic approximation is equivalent to neglecting Hamiltonian terms coupling 
different eigenspaces (see 1 16 1). The evolution inside each eigenspace is, then, generated by 
the following equation: 

= HX t ) - A k (X t )} U k (t), A k (X t ) = ffl fe ^lot (Ai)ri o. (6) 

This equation can be formally solved, and, for a closed loop of the parameters, yields the total 
evolution (notice that by definition 0(X tf ) = 1): 

U k (tf) = T k {t f )V k \ with =Vexp <f A k (X)dX, and T k = £fct , (7) 

where T is an overall dynamical phase factor, and V fe 7 is the celebrated (non-Abeliari) holon- 
omy. In this formula V is the /?af/z-ordering operator, needed because of the non-commutativity 
of the operators A(X) for different values of the parameters. This non-Abelian phases is in 
general very difficult to evaluate, because of the path ordering operation. 



Application to quantum computation 

We would like to mention potential advantages of using geometrical evolution to implement 
quantum gates. First of all, there is no dynamical phase in the evolution. This is because we are 
using degenerate states to encode information so that the dynamical phase is the same for both 
states (and it factors out as it were). Also, all the errors stemming from the dynamical phase 
are automatically eliminated. Secondly, the states being degenerate do not suffer from any bit 
flip errors between the states (like the spontaneous emission). So, the evolution is protected 
against these errors as well. Thirdly, the size of the error depends on the area covered and is 
therefore immune to random noise (at least in the first order) in the driving of the evolution. 
This is because the area is preserved under such a noise as formally proven by DeChiara end 
Palma |4|. Also, by tuning the parameters of the driving field it may be possible to make 
the phase independent of the area to a large extent and make it dependent only on a singular 
topological feature - such as in the Aharonov-Bohm effect where the flux can be confined to a 
small area - and this would then make the phase resistant under very general errors. 
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Figure 1: (a) The geometric phase for a single qubit is proportional to the solid angle Q. (b)The 
four level system that can be used for non-Abelian quantum computation to encode one qubit of 
information in two degenerate levels. The method is detailed in the text. 



So, in order to see how this works in practice we take an atomic system as our model 
implementing the non-Abelian evolution. We'll see that quantum computation can easily be 
implemented in this way. The question, of course, is the one about the ultimate benefits of 
this implementation. Although there are some obvious benefits, as listed above, there are also 
some serious shortcomings, and so the jury is still out on this issue. 



Example 

Let's look at the following 4 level system analyzed by Unanyan, Shore and Bergmann |28|. 
They considered a four level system with three degenerate levels 1,3,4 and one level 2 with 
a different energy as in Fig £>. This system stores one bit of information in the levels 1 and 2 
(hence there is double the redundancy in the encoding of information). We have the following 
Hamiltonian 



H(t) = 
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where P, Q, S are arbitrary functions of time. It is not difficult to find eigenvalues and eigen- 
vectors of this matrix (exercise!). There are two degenerate eigenvectors (with the correspond- 
ing zero eigenvalue for all times) which will be implementing our qubit and they are 

3>i(t) = (cos#f, 0, — sin# t 0) and $2(£) = (sin</) t sin(9 t , 0, sm<pt cos£?t, — cos^t) 



where tan0 t = P(t)/Q(t) and tan</> t = Q(t)/y/P(t) 2 + Q(t) 2 . In the adiabatic limit, we 
can restrict ourself to these states only. Although, in general, the Dyson equation is difficult to 
solve, in this special example we can write down a closed form expression 1281 . The unitary 
matrix representing the geometrical evolution of the degenerate states is 



B(th) = 



cos rjt sin r/t 
— sin rjt cos r\ t 



(8) 



where rjt = /J sint/v^dr. This therefore allows us to calculate the non-Abelian phase for 
any closed path in the parametric space. After some time we suppose that the parameters 
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return to their original value. So, at the end of the interaction we have the matrix B(j]f) 
where rit = £ 9 (SdP — PdS), which can be evaluated using Stokes' 

theorem (the phase will in general depend on the path, as explained before). So, we can have 
a non-Abelian phase implementing a Hadamard gate. With two systems of this type (mutually 
interacting) we can implement a controlled-Not gate and therefore (at least in principle) have 
a universal quantum computer (see [29 1). 
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